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Abstract 

Within the RSJ model, we performed a theoretical analysis of phase-locking 
in elementary strongly coupled Josephson junction cells. For this purpose, we 
developed a systematic method allowing the investigation of phase-locking in 
cells with small but non-vanishing loop inductance. The voltages across the 
junctions are found to be locked with very small phase difference for almost 
all values of external flux. However, the general behavior of phase-locking 
is found to be just contrary to that according to weak coupling. In case 
of strong coupling there is nearly no influence of external magnetic flux on 
the phases, but the locking-frequency becomes flux-dependent. The influence 
of parameter splitting is considered as well as the effect of small capacitive 
shunting of the junctions. Strongly coupled cells show synchronization even 
for large parameter splitting. Finally, a study of the behavior under external 
microwave radiation shows that the frequency locking-range becomes strongly 
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flux-dependent, whereas the locking frequency itself turns out to be flux- 
independent. 
74.50 
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I. INTRODUCTION 



Josephson junction arrays are considered as candidates for microwave oscillators with 
possible applications in the field of satellite communication systems, astronomical observa- 
tions, construction of supercomputer chips and spectroscopyEHI. They are potentially well 
tunable over a relatively wide frequency range while radiating on a narrow linewidth, and 
they can deliver large output power, at least in comparison with a single element. Lin- 
ear arrays of Josephson junctions have been subject to theoretical investigation for more 
than twenty yearJlli'lll^Hli. During last years, there have been some promising experimental 
result&0. Since the beginning of the nineties, there has been a growing interest in 2D 
array^llM. Up to now, the radiation power of 2D arrays was found to be much smaller 
than that of ID arrays (0,1 fiW maximumEl, in comparison to 50 fiW in ID array^l). This 
may be a consequence of technological problems as well as of general properties of 2D arrays. 

Because of this fact, there is some renewed interest in the general mechanisms of syn- 
chronization of coupled Josephson junctions. Most of the early adiabatic investigations 
were based on weak (preferably inductive) coupling of the elements0'&^, which is surely 
fulfilled for relatively large circuits with total inductances ~ 10 pH. However, present day 
technology allows the preparation of very compact arrays with circuit dimensions around 1 
fim and smaller, having inductances smaller than 3pH. In this case, the adiabatic methods 
developed earlier fail. On the other hand, neglecting inductances at allil seems to be a too 
rough approximation. 

In this paper we develop a method for dealing with Josephson junctions coupled via a 
small inductive shunt (SQUID-type coupling) in a systematic way. In Sec. II, we describe the 
circuits handled by this method, derive the general equations, review some results obtained 
by conventional methods and discuss their inapplicability to our problem. In Sec. Ill, we 
give a review of our analytical method and present some results for the simplest case of non- 
hysteretical, identical junctions, which are compared with numerical simulations. Sec. IV is 
devoted to non-identical junctions and Sec. V to the infiuence of an additional small junction 
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capacitance. In Sec. VI we consider the synchronization of strongly coupled SQUID cells 
by external microwave radiation. A summary as well as some speculations about possible 
applications of our results are given in Sec. VII. 

Contrary to most theoretical investigations which are mainly based on computer simula- 
tions during the last years, we concentrate on approximate analytical results. The advantage 
of this approach is usually a better insight into physical mechanisms in combination with a 
broader range of applicability concerning the choice of parameters. The disadvantage is a 
larger amount of mathematical machinery necessary even in relatively simple cases. How- 
soever, we found it quite valuable to use both methods and compare the results. Some 
material described in Sec. Ill was published in a short note earlier0. 

II. A SHORT REVIEW ON WEAKLY COUPLED JOSEPHSON JUNCTIONS 

The strong coupling method described in this paper was developed for the investigation 
of the three similar SQUID-like cells shown in Fig. |l], but the general principles should have 
a much wider range of applicability. These three circuits have a bias current 2Jo, a net loop 
inductance L and a parallel biasing scheme in common. Elaborating the equations of motion 
within the RSJ model, one obtains for the case of identical junctions 



where 0i and 02 are the respective Josephson phases, and iq is the normalized bias current. 



which is supposed to fulfill the condition > 1, here. / marks the normalized loop induc- 
tance. 




(2.1a) 




(2.1b) 



io = lo/Ic {Ic- critical current). 



(2.2) 



/ = 27iIcL/%, 



(2.3) 



and (fi the normalized external flux. 
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(y9 = 27r$/$o {^'- external flux, ^q: flux quantum). 
Dots denote differentiation w.r.t. the scaled time 



(2.4) 



s = {2e/h)IcRN't {Rn- normal resistance). (2.5) 

There exist several investigations of these systems for weak coupling, i.e. / ^ 1 Bill. In 
this case, coupling can be neglected to 0*^ order w.r.t. and both junctions oscillate with 
the Josephson phase of an overcritically biased free contact, 

Co , fCoS-Si/2 



01/2 = arctan 
showing an oscillation frequency 



io + 



— tan 



+ 2 ' (2-6) 



Co = V^o - 1 (2.7) 



and constant phases 6i and 62- To order the lowest harmonics of solutions (|2.6| ) are 
inserted into the r.h.s. of Eqs. |2.14 b. This provides so-called reduced equations for the 



mean values of the phases, averaged over short time scales of the order of Co- Looking for 
solutions of the phase-locking type, 

< 5 > = 0, <S >=< 6i> - <52>, (2.8) 

one finds 

< 6 > r- sin < 6 >= (p. (2.9) 

io[io + Co 



For usual operation regimes with zq ~ 1-5 this results in an approximate linear relation 
between 5 and 0, as indicated in Fig. |^. The normalized voltages 

Vi/2 = (2.10) 

J-CJ^N 



are obtained as 



^1/2 = -- $ (2.11) 

lo + cos(Cos - di/2) 
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This procedure and similar adiabatic methods led to a general understanding of the be- 
havior of weakly coupled Josephson junctions. However, they are not applicable to small 
circuits (of diameter ~ 1/im) with inductances ~ 1 pH. The point is, that for small induc- 
tances the parameter is not longer small, e.g. ^ 3 for L = 1 pH and Ic = 100 /iA. 
In this case the term in parenthesis on the r.h.s. of Eqs. |2.1a| and |2.1b| , resp., dominates the 



remaining terms and it is not possible to derive reduced equations of the type mentioned 
above. 

In order to deal with small inductances which already can be realized technologically 
today, we developed an alternative systematic scheme, which is described in the next section. 

III. A PERTURBATION SCHEME FOR STRONG INDUCTIVE COUPLING 

At first, we found it convenient to introduce new variables 

A = (02-0i)/2 and S = (02 + 0i)/2, (3.1) 

providing the set of equations 

A-FcosSsinA = y((/?-2A), (3.2a) 
E + sin S cos A = zq- (3.2b) 

These equations already indicate that the behavior of S is determined necessarily by the 
bias current io and that of A by the external flux, although coupling makes things slightly 
more complicated. 

We perform a perturbation expansion valid for small /, 

A = Ao + /Ai + 0(/2), (3.3a) 
E = So + ZSi + C(/2). (3.3b) 



This ansatz resembles the " slowly- varying amplitude method" developed several years agoEJcll 
for then and resulting in reduced equations similar to those mentioned above. Expanding 
the sin and cos terms, it is not sufficient that the condition 
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I^Sil < |So| (3.4) 

is fulfilled, but additionally we must demand the much more rigid conditions 

and |/Ai|<l. (3.5) 

Investigation of the final result shows, that this can be met for all (p by choosing a sufficiently 
small /. Only the case ip = tt requires a slightly more involved consideration. 

After inserting the expansion ( p.3a| ,b) into Eqs. |3.24 b one can compare identical powers 
of /. To lowest, i.e. (—1)**, order one finds 

Ao = p/2 (3.6) 

without solving any differential equation. This general feature is valid for higher orders of 
/, too: The variables S„ must be determined by solving a first-order differential equation 
(only to lowest (O*'^) order being non-linear), whereas the A's can be calculated algebraically. 
Introducing a once more rescaled time, 

s = s-cos{p/2), '=^, (3.7) 

ds 

the equation for Sq (O*'^ order w.r.t. /) becomes very similar to that of an autonomous 
junction, 

So + sinSo = io, "^o = '^o/cos(v3/2). (3.8) 



Thus, we can immediately write down the solution 

' Co 
io + cos(v3/2) 



2arctan ( — ^;%--^tan(^)^ + ^, (3.9) 



with 

Co = v/^§-cos2(^/2) (3.10) 

where we have imposed the initial condition 

So(s = 0) = 7r/2. (3.11) 
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Because the A„ are determined algebraically there is no freedom to specify separate initial 
conditions for them. To ensure the validity of the perturbation expansion, it is important 
to specify all higher S„ according to 



= 0) = for n>l. 

With Aq and Sq given Ai can be determined algebraically, 

1 / Cosin(Cos) 



(3.12) 



Ai 



sin((/9/2) 



2 + cos((/9/2) cos(Cqs) 
This expression automatically satisfies the initial condition 

Ai(0) = 



(3.13) 



(3.14) 



required for the validity of the perturbation expansion. For Si one obtains the inhomoge- 
neous linear differential equation 



Si + Si cos So cos(v9/2) — Ai sin Sq sin(y9/2) = 0, 



(3.15) 



admitting the solution 

tan2(^/2) 



Si 



^ I /o^n 7 N , 72, ^o + cos(y/2)cosCog ^ 
2ocos(v?/2)(l - cosCo-s) + CqI^^ ■ 7 — TTTn (3.16) 



2(io + cos((y9/2) cosCo'S) V" — ~,u"y ' 2o + cos(v9/2) 

where we already exploited the initial condition, Eq. p.l2| . From Eqs. pl9| , p.l3| , and 



3.16 one derives the normalized voltages 



Vl/2 



+ / 



Cosin(v9/2) 



%Q + cos((y9/2) cosCoS 2(^0 + cos(v9/2) cosCo-s)" 



(3.17) 



■ I in\ ■ 7 { ■ ^ t /o\ > Co 1 zq + cos y?/2 cosCpg 

sm^/2 sm (qS «o + cos y?/2 + — — In — -— 

\ cos(v9/2) lo + cos(v9/2) 

=FCo^cos(9?/2) + io cos Cos 

In this way we derived an analytic expression for the voltage on junction 1 and 2 resp., 
valid to first order w.r.t. /. It has to be compared with the solution ( |2.11D for weakly coupled 
elements. Both have in common, that there always exists a solution showing phase-locking 
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for all values of the external field. However, contrary to the case of weak coupling, where the 
phase shift scales approximately with external flux, we will shortly see that it is negligibly 
small for all values of the external field, except for ^ tt. On the other hand, the frequency, 
being flux-independent for weakly coupled elements, becomes flux-dependent according to 
Eqs. |3.17| ,10, detuning the cell this way. 

So far, we could not derive an equation of motion (a kind of reduced equation) for the 
phase shift of strongly coupled elements directly from the basic equations as it is possible 
for weakly coupled elements. Therefore we must adopt an appropriate alternative definition 
of phase difference. We define phase shift as the difference of the mean value crossings of 
the lowest harmonics of Vi and V2 according to ( p.l7| ). 



To proceed, we must evaluate the lowest Fourier coefficients of Eqs. |3.17 



o; — — 

vi/2 = - + ai/2 cosCqS + bi/2 sinCoS. (3.18) 

After some calculation including a Taylor expansion of the logarithm one obtains 

a = 2Co, (3.19) 

ai/2 = -^(-2 cos(^/2) T Ko sin(¥^/2)), (3.20) 
«o + Co 

, ,^,i^im^^L^s^^^^_\. (3.21) 

N + Co) V '0 "^^oi^o + Co)) J 

Note that the Fourier coefficient b being proportional to / is small compared to ai/2; the 

coefficient ai/2 itself is dominated by its first term, except for (p = tc. Thus, with the possible 

exception of the vicinity of this value, both voltages are in phase independent of the external 

flux. 

Adopting the definition given before, one can derive a formula for the phase shift 5 as a 
function of the Fourier coefficients, 

5 = arccos ( -_ii^±|=_ ) , (0 < ^ < tt). (3.22) 

Deriving this formula, no assumption has been made about the structure or order of magni- 
tude of the Fourier coefficients. Extension to n < ip < 27t needs a more subtle investigation 
of the solution itself to get the correct branch of the arccos; in this case one obtains 

9 



6 = sgn((/9 — tt) 



TT — arccos 



aia2 + 6^ 



b^{al + al + b^)+alal 

(0 < < 27r) 



vr, (3.23) 



Fig. 1^ shows the phase shift between f i and V2 as a function of external flux (p. 

Our analytical approach was accompanied by numerical investigations exploiting the 
Personal Superconducting Circuit ANalyzer program PSCANi§ (Fig. H). Comparison 
of Figs. ^ and | shows that even for / = 1 where the analytical approximation is not 
longer valid results are quite similar to those of the numerical simulation. Both figures show 
that already in this case the behavior is qualitatively different from that of weakly coupled 
elements. Thus, in the intermediate regime / ~ 1 the strong coupling scheme provides a 
better approximation to the actual behavior of the SQUID cell than the weak coupling 
results do. 

There are two limiting cases of special interest. For vanishing external flux the elements 
behave in the same way as one junction; the phase shift between voltages vanishes and the 
elements oscillate with voltages 

Vi = V2= . , ^° . . (3.24) 
zo + COS Cos 

This behavior is plausible, because in this case there is no current flowing through the shunt 
and both elements act essentially as one element. For ip = tt one obtains 

Hq I 

vi = io — — sinios + -sm2iQS, (3.25) 

Hq I 

V2 = io + sin ios + - sin 2ios. (3.26) 

This clearly indicates that the phase difference S between the voltages according to our 
definition is equal to vr. Within the narrow transition region between the regimes 5 ~ and 
6 ^ IT oscillations become highly non-harmonic (cf. the voltages in Figs. ^ and ^). This is 
caused by higher harmonics becoming dominant in comparison to the lowest harmonic with 
frequency Co- this case the definition of phase difference based on the lowest harmonic 
of the Fourier expansion (|3.18|) as well in the analytical calculation as in the numerical 
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procedure becomes less reliable. The jump observed in the numerical curve (Fig. ^) may be 
assigned to this fact. 

To derive a more simple rule of thumb Eq. |3.22| can be rewritten as 



5{ip) = arctan(a2/6) — arctan(ai/&), < ip < n. (3.27) 

with ai/2 and b according to ( p.20| , p.2lD . We consider this formula for very small, but finite 
/. In this case, the second arctan approaches the value —7t/2. The first one changes its sign 
at the flux (f*, 

cos{ip*/2) - ^Co sm{ip*/2) = 0. (3.28) 

Exploiting this formula while neglecting higher orders in I, one obtains 

ip* ^71 - iol. (3.29) 

This provides a simple approximation for the phase shift of the cells under investigation, 

5 ^ 7r0((/? - (/?*) (y?<7r). (3.30) 

Fig. 1^ shows that for sufficiently small / the solution is indeed perfectly approximated by a 
Heaviside step function. This approximation might be useful considering more complicated 
arrays. 

Another quantity of interest is the I-V characteristics of the cells under investigation. 
From (|3l8| ) one easily obtains 



Vy2 = ^il - cos2((^/2). (3.31) 

This reproduces a well-known resultH: The l-V characteristics of a small-inductance SQUID 
has a hyperbolic shape, the vertex being dependent on the external flux. 

IV. PARAMETER-SPLITTING IN STRONGLY COUPLED CELLS 



Real junctions never have identical parameters. The response to parameter differences 
becomes particularly important in large arrays, where one usually cannot avoid a parameter- 
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splitting of the order of one p.c, at least. In this section we consider junctions having 
different critical currents as well as normal resistances, 

-^Ci 7^ Ic2, 7^ Rn2, (4.1) 

with the subsidiary condition 

IciRni = IC2RN2, (4.2) 

which is usually realized as a consequence of the technological process with a good accuracy. 
Introducing the mean critical current 

Ic = lilc, + IC2) (4.3) 

and the parameter-splitting 

^ = 7^^. (4-4) 

one derives the following RSJ-model equations for the cells shown in Fig. |1|: 

01 + sin 01 = - (01 - 02 + ¥?), (4.5a) 

02 + sin02 = Y^ + I^Y^ (01-02 + <^). (4.5b) 



As before, it is advantageous to introduce new variables A and S according to Eq. |3.1| . 
Eqs. p. 24 b are then modified to 



A + cos E sin A = zo + J^^^^iv - 2A), (4.6a) 

S + sin S cos A = ^ - J^^^^iv - 2A). (4.6b) 

Some effects are already qualitatively displayed by this couple of equations, (i) To first order 
in ■}) there is a correction of the magnetic flux ~ —ioW. (ii) There is a correction of the 
bias current i^lip — 2 A) /I being of first order, too. It includes an additional coupling via 
A. Eqs. ( [4.64 b) indicate that for weak coupling {I ^ 1) the additional magnetic 

flux dominates, an effect which has already been observed (cf. Eq. 13.30b ini). However, 
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for strong coupling (/ <^ 1,^9 <^ 1) this term is of second order only. It turns out that the 
difference (y? — 2 A is of the order of /, so the addition to the bias current is of first order in 
and dominates. 

First of all, we are interested in the maximum parameter-splitting which is possible 
without destroying synchronization. For this purpose, the splitting-parameter should not 
be considered small from the beginning. As before, we perform a perturbation expansion 
w.r.t. / according to ( p.3a| , ^^3b| ). To lowest, i.e. —1^*, order we again obtain Aq = To 



order, we get the system of equations 



So + sinSoCos(v3/2) = + ^^^Ai, (4.7a) 



■}) 2 

cos Eo sin Ao = - j^^^o - Y^^^' ^^'^^^ 

Again, Sq has to be determined by solving a differential equation, whereas Ai is calculated 
algebraically. The new feature is an additional coupling between both variables caused by 
the last term on the r.h.s. of Eq. |4.74 Combining both equations, one obtains 



Eo + sin So cos(v9/2) + cos So sin(v9/2) = ^ (4.8) 



In comparison to (|3.8|) this equation shows an additional non-linearity due to the parameter- 
splitting. Moreover, the current is divided by (1 — -i?^), indicating that a splitting of critical 
current can lower the current necessary for the onset of oscillations. 

Eq. ^]8| can be handled exactly. There are four different types of solution^. Only one 
of them shows the continuous transition to the case = and the corresponding oscillating 
voltage: It is realized for 

'° > sin2(^/2) + cos2(¥;/2). (4.9) 



1 

Further estimation gives the bound for oscillations to occur, 

il>l-^^. (4.10) 

For = this reproduces a well-known fact. In addition, it proves the conjecture above 
that growing parameter-splitting leads to an effective lowering of the critical current. 
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With this condition fulfilled we could evaluate Sq and from ( |4.7b|) Ai. However, although 



Si has to be determined from a first-order linear differential equation, the resulting integrals 
are rather intricate. Thus, we performed a perturbative treatment of the system ( ^.6a| ,b) 
not only w.r.t. /, but w.r.t. i), too. This is more delicate, of course, because there are 
two parameters involved. To discuss the /- and '(9-dependence independently, its not wise to 
specify the ratio l/'d from the beginning. We only suppose / <^ 1 and ^ <^ 1, leaving the 
ratio -^/l unspecified. To first order, we write down the expansion 

A = Ao + /A10 + ^9Aoi, (4.11a) 
S = So + /Sio + 79Soi. (4.11b) 

Inserting into ( [4.6a| ,b) and comparing equal orders one obtains the set of equations 

necessary to evaluate the A's and S's. For Aq, Sq, Aio, and Sio one obtains similar equations 
as before. Furthermore, one observes Aqi = 0; thus, no additional phase shift is caused by 
the parameter-splitting. For Sqi one obtains an equation similar in structure to that for Sio 
(cf. Eq. |37[5| ), 



Soi + ^01 cos So cos(v3/2) + cos Sq sin(y3/2) = 0. (4-12) 
It admits the solution 

Soi = . J sin(^/2). (4.13) 
iQ -|- cos((^/2) cosCqS 

Thus, weak parameter-splitting leads to an additional in-phase contribution 

^(^0 + cos(^/2)) sin Co. ,. (4.14) 

(«o — cos(v9/2) co^QqsY 



to be added to the voltages (|3.17|) . The Fourier coefficients ai/2 are unaffected by this. 



whereas there is an additional contribution to 61/2, 

6f/2 = 2^(.o + cos(^/2))^^^^. (4.15) 

The solution obtained this way proves our earlier conjecture on the dominant contribution 
in the strong-coupling case (Fig. |^). One observes that the phase shift, being slightly raised 
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generally is considerably lowered for ip = t[. To lowest order of our analytic approximation 
(valid for strong coupling and weak parameter splitting) there is no indication of a shift of 
the peak caused by the parameter splitting. This is confirmed by comparison with numerical 
simulation, as long as parameter splitting is sufficiently small (t? ~ 0.2). For larger the 
numerical result (Fig. |^) gives a first hint to the peak shift. 

V. CAPACITIVELY SHUNTED JUNCTIONS 

The infiuence of the displacement current fiowing through the junctions was neglected 
up to now. This is justified, as long as the McCumber parameterSS 

(3 = '^IcR%C (5.1) 

is negligible. 

In this section, we will investigate the infiuence of /3 7^ 0. The displacement current 
adds a second-derivative term to the RSJ model equations (sometimes called RCSJ model 
equations then), 

+ 01 + sin 01 = io- l~^{4>i - 4>2 + V^), (5.2a) 

P<p2 + 02 + sin 02 = io + r\(pi - 02 + <^). (5.2b) 

In general, the second derivative may change the character of the differential equations 
completely; for instance, it is well-known that there appear new types of solutions showing 
chaotic behavior0'ii. Here, we will restrict our treatment to small (3 (/3 <^ 1) guaranteeing 
a continuous transition to the former solution for /9 = 0. 

Again, it is recommended to combine Eqs. ^.24 b obtaining 



/3A + A + cosSsinA = y((^-2A), (5.3a) 
+ S + sin S cos A = io- (5.3b) 

One clearly sees, that both equations are affected by the additional /?-terms. Again, both / 
and P are supposed to be small parameters justifying the expansion 
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A = Ao + /Aio + /9Aoi, (5.4a) 
E = So + /Sio + /3Soi. (5.4b) 

The resulting equations for Aq, So,Aio, and Sio, resp., are essentially the same as before. 
For Aoi one readily recovers 

Aoi = 0. (5.5) 

The only new equation concerns Sqi, 

Soi + Soi cos So cos (v9/2) = — So, (5.6a) 

where we have already exploited some of the previous results. Again, this is an inhomoge- 
neous linear differential equation, with the inhomogeneity being determined by the already 
well-known So. The solution, obeying the correct boundary condition (Soi(s = 0) = 0), is 

Soi = /° , - In ^o + cos(^/2)cosCoS _ 

io — cos{(p/2) cosCqS zo + cos(v9/2) 

One obtains a logarithmic structure similar to that observed earlier in formula ( p.l6|) . The 
term ( p. 7] ) provides to both voltages a contribution 



. ^ Co cos(y)/2)sinCog A + cos(v9/2) cosCo-s _ 
(io + cos(v9/2) cosCos)^ V io + cos{(p/2) ^ 



(5.8) 

Because the logarithmic structure is already present in ( |3.17| ), it is not hard to evaluate the 
corresponding capacitive contribution to be added to the Fourier coefficients bi/2 according 

6f,. = ''"=-2^CS^^^^(l + ^-|f^). (6.9) 

(^0 + Co) V ^0 ^^o{^o + Co) / 



The phase difference obtained by inserting (|3.21|) and ( |5.9|) together with the unchanged 



coefficients ( |3.20| ) into ( |3.23| ) is shown in Fig. Results of a numerical calculation per- 



formed in parallel are given in Fig. Qd. The general tendency is that a non-vanishing 
capacitance (/5 ~ 1) slightly enhances the phase shift without qualitatively changing the 
general behavior. For (3 > 0.5 the agreement between analytical approximation and numer- 
ical simulation becomes worse, but the same general tendency is still preserved. This is, of 
course, simply a result of the fact that higher orders in /3 are no longer negligible. 
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VI. STRONGLY COUPLED SQUID CELLS UNDER MICROWAVE RADIATION 



There is some interest in the behavior of the SQUID cells under microwave radiation 
from at least two points of view. First of all, the topic is interesting for the construction 
of sensitive microwave detectors. Secondly, knowledge of the behavior under microwave 
radiation is necessary for the study of synchronization in larger arrays, where the long-range 
interaction via external shunts acts similar to an external microwave signal. 

The external microwave signal can be described by an additional ac current, leading to 
the system of equations 

+ sin (f)i = io - - (j)2 + (p) + iujsmujs, (6.1a) 

02 + sin 02 = io + ^""^(01 - 02 + V') + ^o; sincjs. (6.1b) 

As a result, only the equation for the sum variable S is affected and becomes 

S + sinScos A = io + ^ojsincus, (6.2) 



whereas Eq. |3.2a| for the difference variable A remains unchanged. We apply a perturbation 



scheme similar to that used before. From the beginning we will assume / <^ 1, as before, 
justifying the expansion ( p.3a| ,b). Solving for Eq and Si, resp., we must take i^j to be small 
in some intermediate steps, too. 
The first result is 

Ao = ip/2, (6.3) 

as usual. The corresponding equation for Sq, 

So + sin So cos{ip/2) = io + i^j sinus, (6.4) 

is decisive for the behavior of the solution. Introducing the scaled time s according to (p. 7]), 
we obtain 

So + sinSo = io + z<^sina;s, = t^-jt^, = r-^. (6-5) 

cos(v9/2) cos(v9/2) 
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which formally has the same structure as the equation describing an autonomous Josephson 
junction under external irradiationllMi. 

It is well-known that phase-locking of an autonomous junction takes place only if the 
frequency of the external microwave does not deviate too far from the inherent Josephson 



frequency Co = V^o - 1, 

< 1 (autonomous contact). (6.6) 



2C0/ > ^ 

— (^-Co. 



The main new feature in our case is, that the corresponding quantities substituted for io,iuj, 
and u resp. according to 

io io, (6.7) 

iuj iuj, (6.8) 

uj LJ (6.9) 

are dependent on the external magnetic flux. Exploiting the corresponding equation 

< 1, (6.10) 



2Co / ~ ; X 
— (^ - Co) 



one obtains the phase-locking condition 

Co - ^^^^ < - < Co + (6.11) 

An interesting conclusion from Eq. |6.11| is that for external flux cp = n the range of phase- 
locking shrinks to a single point, u = = i^. In this case, for all practical purposes 
phase-locking disappears at all. This is confirmed by Fig. |^, showing the range of phase- 
locking against external flux (f. The reason for this behavior is obvious from examining ( |6.4D : 
For (f = IT the non-linear term vanishes thus removing the non-linearity of the equation at 
all. One more observation is that the center of the locking-range determined by Co becomes 
fiux-dependent, too. 

In case of phase-locking, phase shift between the external microwave and the circuit (the 
last being characterized by the sum variable S) can be deduced from comparison with the 
autonomous contact yielding 



^o = arcsm ^^°^^^ (6.12) 
luj cos(v9/2) 

To verify our result in an independent way, we performed a numerical operation range 
analysis automatically integrating the corresponding differential equations and checking, 
whether the results lie within a certain bound. The output of this analysis indicated by 
diamonds in Fig. |^ is in excellent agreement with the analytical results. In view of the 
experimental setup the figure should be interpreted as follows: For a fixed bias current and 
frequency of the external microwave radiation there is a limited range of flux indicated in 
Fig. ^ within which synchronization occurs. Within this range the whole cell oscillates 
with the microwave frequency, uo. For small microwave intensities this puts rather severe 
conditions on the external flux, as is observed by comparing Figs. b and c. 

In case of phase- locking, we obtain to 0*^ order of perturbation theory 

Eo = 2arctanf tan f — - + - (6.13) 

" V^o + cos(v9/2) V 2 2 2 ^ ^ 

with the lowest-order phase shift 5q according to ( |6.12|) . 



Within the next perturbative order, Ai is determined algebraically as before. The so- 
lution is identical to Eq. |3.13| , the only difference being, that within the time-dependent 
arguments one has to substitute 

Co-s ojs — 5q. (6.14) 

For El we obtain the equation 

Si + Si cos So cos(v9/2) — Ai sin Sq sin(y9/2) = 0. (6.15) 



Substituting Sq we exploit ( |6.4| ) neglecting the external current, bearing in mind that 
is small and Si is already of first order. Within this approximation, the solution has the 
same general structure as ( |3.16|) , where we again have to substitute ( |6.14[ ). As a result, in 



addition to the fiux-dependent phase shift between the SQUID circuit voltage oscillations 
and the external microwave signal we obtain the same (mostly negligible) phase splitting 
between the junction voltages than without radiation. 
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To summarize, solution (|3.17|) is reproduced with the only substitution ( |6.14| ). This has 
several consequences: (i) The frequency of the oscillations is determined by the microwave 
frequency only and turns out to be independent of external flux within the locking-range, 
(ii) If external flux is present in addition to the microwave radiation, this flux will limit the 
range of phase-locking in general and destroy synchronization at all in case of (f = tt. (iii) 
The relative phase of both junctions is not influenced by the external radiation up to flrst 
order in perturbation theory w.r.t. /. (iv) There is an additional shift of both phases relative 
to the external radiation according to Eq. |6.12| , which is controlled by external flux as well. 



VII. SUMMARY 

We investigated the synchronization behavior of three similar 2-junction SQUID cells 
with strong inductive coupling. For this purpose we developed a perturbation scheme ap- 
propriate for small but non-vanishing inductances. The perturbation ansatz itself is in a 
certain sense similar to the " slowly- varying amplitude method" developed several years ago. 
However, application to strong coupling completely changes the character of the expansion. 
Generally, the procedure is more involved than in the case of weak coupling. Therefore 
we were not able to derive an explicit equation of motion for the phase difference between 
voltages. In view of this fact, we determined voltage phase shift from the lowest Fourier 
coefficients of the voltages. 

For identical junctions without capacitive shunting we find that for every value of the 
external flux a phase-locking between oscillating voltages takes place like in weakly coupled 
elements. However, contrary to the case of weak coupling, the phase shift is negligibly small 
for almost all values of external flux, with the exception of the vicinity of if = it. On the 
other hand, the frequency, being flux-independent for weakly coupled elements, becomes 
strongly flux-dependent in strongly coupled elements and, consequently, the corresponding 
I-V characteristics too. The results obtained are compared with numerical calculations. 
Generally, a good agreement is observed. Especially it is found, that the strong-coupling 
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approach provides a good approximation not only for very small inductances I, but is better 
suited to describe the intermediate range I ~ 1 than the ordinary weak-coupling approach. 

If both junctions are not identical the influence of parameter splitting is found to be 
qualitatively different for weak and for strong coupling. For weak coupling, parameter 
splitting mainly leads to an additional contribution to the external flux; as a result, the 
whole phase-flux-dependence is shifted by some value. In case of strong coupling, this effect 
can be neglected and the leading contribution is a correction of the bias current. This 
correction acts in favor of synchronization and lowers the phase shift present in a small 
range around ip — tt. 

In case of identical junctions having a small, but non-vanishing capacitance the main 
result is a slight enhancement of the phase shift, although the qualitative picture is not 
changed, at least if /3 < 1. For /3 > 0.5 the agreement of analytical results and numerical 
simulation becomes less convincing, obviously showing the limitations of applicability of the 
analytical perturbative method. We should mention that in the case /3 7^ as well as for 
parameter splitting two independent expansion parameters must be considered small. 

Finally, we investigated the behavior of the cells under external microwave radiation. In 
this case, we observed a limited locking range similar to that of an autonomous Josephson 
junction under external radiation. However, for a strongly coupled cell the synchronization 
range is flux-dependent and shrinks to zero for (f — ir. In addition, the width of the 
synchronization range depends on the amplitude of the external radiation. 

Contrary to our results for a freely oscillating cell, external radiation synchronizes the 
cell in such a way that the oscillation frequency becomes flux-independent within the flux- 
dependent locking range. However, within the synchronization range, an additional phase 
shift between external radiation and internal oscillations takes place as well as a shift of the 
whole synchronization range. 

From our study one can draw the general conclusion, that two strongly inductive coupled 
junctions behave like a free junction if no flux is present within the cell. External flux 
tends to shift the phases between the voltage oscillations of the two junctions, but for most 
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practical applications this splitting is negligibly small. More serious is the fact that already 
in a cell of identical junctions the oscillation frequency itself becomes field-dependent. This 
detuning of the cell has several consequences for the construction of larger arrays. First 
of all, external fluxes must be shielded, preferably by an external superconducting ground 
plane. Secondly, additional fluxes are produced by the array itself, which might seriously 
disturb the synchronization^. 

A possible way to circumvent this problem and to obtain phase-locking could be to include 
an external long-range interaction via an additional shunt. The methods developed in Sec. 



V| could be helpful in such an investigation. For instance, according to our observations, 
an external shunt with sufficiently strong coupling-strength to make the synchronization 
frequency flux-independent may play a crucial role for obtaining phase-locking in large two- 
dimensional arrays. 
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FIGURES 

FIG. 1. Three SQUID Cells which can be described with the strong coupling method described 
in this paper. 

FIG. 2. Mean voltage phase shift 6 against normalized external flux ip from analytical approx- 
imation for weak inductive coupling (iq = 1.5). 

FIG. 3. Phase shift 6 against normalized external flux ip for strong inductive coupling I = 0.1 
and medium inductive coupling I = 1.0 obtained from analytical approximation (|3.23| ) (zq = 1.5). 

FIG. 4. Like Fig. results obtained from numerical simulation. (A tiny shunt capacitance 
P = 0.01 has to be added here.) 

FIG. 5. Time-dependence of the voltage for two different values of the external flux, a) (p = 7r/2 
h) if = IT, Parameters: = 1-5, I = 0.1. 

FIG. 6. Phase shift 6 against normalized external flux (/? for extremely strong coupling 
(0.001 < / < 0.1), obtained from analytical approximation ( 3.23| ) (iq = 1.5). 



FIG. 7. The influence of parameter splitting on the phase shift 5 against normalized external 
flux (f for strong coupling obtained from a) analytical approximation ( |3.23D , b) numerical simulation 
(io = 1.5, / = 0.1, 2 < 99 < 4). 

FIG. 8. The influence of a non-vanishing capacitive shunt of the Josephson junctions on the 
phase difference for iq = 1.5, / = 0.1, f3 = 0.2, 0.5, 1.0. a) analytical approximation, b) numerical 
simulation (io = 1.5, / = 0.1, 2 < ip < 4). 

FIG. 9. Synchronization range for a strongly coupled SQUID cell under external microwave 
radiation with frequency u> for a)ia, = 0.1, b)*^^ = 0.2, c)iuj = 0.4 (bias current zq = 1.5 ). The 
diamonds indicate results from numerical operation range analysis. 
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